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The study is concerned with the interference of friction in constant temperature heating inside the diverging part of a converging- 
diverging nozzle working at design condition where the velocity changes to supersonic flow regime from a Mach no. of 1.0. Set of 
gas-dynamic equations are completed using numerical approximation of the dimensionless derivatives. A finite difference technique 
with 20% under — relaxation factor is applied. Most of these equations are non-linear along from the throat till the nozzle exit with a 
wide range of frictional forces from 0.001 to 0.003 average frictional factor ranges with 0.00025 steps. The results are distributed 
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between 3.0 and 4.0 area ratio of non-linear inside nozzle curvature. The multi mixing of all variables help to give results different 


and more accurate values than isentropic or with isothermal or frictional states or with double mixing of any of them to be more 
close to experimental ones. These all together have been considered in our model which has employed all variables which will assist 
in the optimum design for real type of nozzles. 


Keywords: frictional, constant temperature heating, Supersonic air flow, diverging part of nozzle. 


NOMENCLATURE 

A Cross-sectional-area in m? 

Ae Exit cross-sectional area in m2 
Ne Critical cross-sectional area in m? 
a,b &c Area shape constants 

DR Mass density Ratio 

f°) Mass density in kg-m>3 

oe Critical mass density in kg-m-? 
Vv Flow velocity in m-s"! 

v Critical velocityin m-s"! 

VR Velocity ratio (v/ v’) 

ER Area-ratio (A/ A’) 


ERE Exit area-ratio (Ae/ A’) 


f Average friction factor 
F Frictional force in N 
k Specific heat ratio 
L Real length of duct in m 
m° Mass-flow-rate in kg-s™" 
M Mach number 
P Pressure in kpa 
e Critical pressure in kpa 
Po Stagnation pressure in kpa 
Ps Critical stagnation pressure in kpa 
PR Pressure ratio (P/ P’) 


PZR Stagnation pressure ratio (Po/P>) 
Under relaxation factor 


Ga) ou 


Gas-constant in J-kg"'-K"" 
Entropy property in J-kg"'-K"! 


* 


Critical Entropy inJ-kg™!-K"! 


* 


nn 


Entropy change in J-kg"'-K"" 


~ 
v) 


Entropy change ratio 
Temperature in K 


* 


al Way Way ny 2) 


Critical temperature in K 


| 
fo} 


Stagnation temperature in K 


* 


Critical stagnation temperature in K 
TR Temperature ratio (T/ T’) 
TZR Stagnation temperature ratio (To/TG) 


° 


x Distance along nozzle in m 
X Dimensionless distance (x/L) 
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1. INTRODUCTION 


Most of related studies about engines of aircrafts or space-shuttles or different kinds of rockets are concentrated on flows in 


supersonic regime of properties. The assembly of most variables take place will close the studies to a more reality, as long as the 
prediction is coming to exact points. This is very important in modern design condition of nozzles operation. Reference [15] has 
described the experimental study of a supersonic interval flow through a converging —diverging nozzle. The Mach no. distribution 
obtained using traversing of a pilot-Pitot -Tube. The used converging-diverging nozzle has an area-ratio of 2.867. Even the scope of 
their study was to investigate the improvement of flow using types of riblets, but the smooth surface tests can be compared with 
present study. Therefore their study lacks to theoretical analysis, mainly in the supersonic part between throat and nozzle exit.Since 
that study employed a draft wind tunnel, however can be stated as a constant temperature (isothermal) type tests. Also smooth 
surface tests [15] assumed a coefficient of friction of 0.0028 in their tests. They examined the distribution of both the static and 
dynamic pressures. 

Reference [4], has verified and validated of the quasi-one dimensional pressure based on finite volume algorithm, implemented 
in Generalized Fluid System Simulation Program (GFSSP), for predicting compressible flow with friction, heat transfer and area 
change. Their converging-diverging nozzle configuration area ratio was 4.0 with an internal wall average friction factor of 0.002. They 
concentrated on the chocking condition of nozzle mode of operation, i.e. subsonic flow regime within the diverging part of nozzle 
after chocking (M= 1.0) at around throat due to all friction, heat transfer and area. So they left the supersonic mode of operation 
(design condition). 


Area Ratio 
3.0 4.0 


Figure 1 Schematic of a converging-diverging nozzle. 


As well known the inside wall of variable area duct will subject to friction in spite the lower distance along the duct in addition to 
isothermal or heat-transfer induced together. The one of preferred methods to analyze the compressible fluid flow regimes are 
numerical methods. The real power of the numerical approach is, however, felt when problems for which it is difficult to obtain 
analytical solutions are considered. The methods used to analyze one variable separately. Studies on area-variation or friction or 
isothermal are proceeded before in constant area ducts or a mix between two of them. In the present study the analysis with full 
effect of the three types of flows will employed here using a finite-difference numerical model after experiencing it to be more 
applicable with such complex algorithm forms. At the same time a comparison is made with two flow types solved with the same 
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numerical method. Not only this but also examine different real frictional forces in addition to real two area ratio configurations and 


analyze all flow properties in addition to Entropy change, as shown in figure (1). The numerical model deals with the supersonic flow 


regime in the design operation mode condition. The analysis is difficult to employ using normal analytical solutions. 


2. AREA VARIATION FLOW MODEL 


The conservation of mass between any static and critical states in an isentropic area variation model are [1], [7], [11], [16] and [17], 


m°=p-v-A=p*-v*-A* 


Conservation of momentum between same states before is, 
F+P*.A*—-P.A=m(v—v*) 


The fluid force is, 


* 


F=(P°+——)(a-a 


Equation (3) becomes, 


* 


(P+ ) (A= AY) + PLA PLA = p-v?-A- pt ev?’ 


i ASC hae Peep hoes 
5P*. 5P. 5P. 5P* pv p*:v 

1 1 

5A Pie A(R =P) = peed pty ane 


1 1 
sPAt AY) —SP(A+A) = piv? A— pt ev? AY 


* 


seat ayy +( ) ne ar =5P(A+a‘) + (2) 2A 
m1 i ane RT) * 


Inserting Mach no. definition, 


v=M.vkKRT 
A+ A* A+ A* 
P*/ +k.A’) = P( +M?.K.A) 


The non-dimensional pressure ratio will be, 


i(i+a)+k \w 


Be ~(1+2) +k.M? A 


Substitute for (p*/ p) from equation (2), 
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T P.v.A 
TY Pee 
Again using the Mach no. definition, 
T  P.A.MVKR.T 
T* ps. AtVvkK.R.T* 
T  (P.A.M)? 
T* (P*.A*)2 


Substitute for (P/P*) from equation (12), 


r_MmG(tte)+¥) 


T fi At 2\2 
(-(1+4)+kM?) 
The non-dimensional velocity ratio can be obtained from equation (15), 


vi T (=) 
ve TY PA 


Substitute for non-dimensional temperature and pressure ratios from equations (12) and (18), respectively, 
2[1 A 
A eis Gi+n)+*) 
ve oi AY 2 
*(1+2)+kM 


Using equations (2) and (20) to find the non-dimensional density ratio, 


p A (G(1+5)+kM?) 


* awe (2(444 
py a.m? (2(1 +4) +k) 
Then applying the adiabatic law to find the non-dimensional stagnation temperature ratio, 


Ty T(4 +i (k- 1)M?) 


T T(1+&-n) 


Again using equation (18), 


TM? G(a +=) +k) (1 +2(k- 1)M?) 


To (:(1+%)+km2)'(1+20- ») 


The stagnation pressure ratio can be obtained from, 


The non-dimensional Entropy change will be, 
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7 Tyiei 
. = =In (=) (26) 
ae 
Substitute for (P/P*) and (T/T’) from equations (12) and (18), respectively, 
ke 
ss GOP (GOs +) 
R (E(1+4) +k) " *(1+") +k.M? a 
- Ae 2 A 


Note that the stagnation properties and Entropy are constant for isentropic flow only varied for other types of flow as will be shown. 


3. DIMENSIONLESS FORMULATION 


Continuity equation can be repeated using dimensionless values between initial state (1) and (i)th state as follows [11], 


DR(1). VR(1). ER(1) 


DRO = Tay ER@ =-(28) 
Pressure ratio in dimensional form is written as follows, 
m-( *(1+ER) +k Je a 
7 ; an is 
=(1+=:) +k.M? 
The pressure ratio at state (1) is, 
By 
=(1+ER(1))+k 1 
PR(1) = (ee (=a) .»- (30) 
Z os ER(1 
*(14+ aa) +k M(1)2 (1) 
In the same manner the pressure ratio for (i)th state is, 
1 
=(1+ ER(i)) +k 1 
PR(i) = f cole ( ) (31) 
: ER 
5(1+ gap) + M@?] ERO 
Dividing equation (31) by equation (30), 
1 ; 1 1 
PR(i) = PR(1) (i + ERG) + ) 4 Pe) te sar (ERG) (a0) 
_ t 1 1 : : ee 
(1+ ERQ)) +k/\F(1+ 55) +k MG? / ERO 
Velocity ratio in dimensional form is written as follows, 
dPR 
kVRdVR=—+ (33) 
Then make a backward property difference, 
1 1 
VR(i)(VR(i — 1)) = — (= a) () (PR(i) — PR(i- 1)) (34) 
VR(i) = VRGi— 1) —- (Groweo) (PR(i) — PR(i— 1)) (35) 


Samheri, 
Frictional and Constant Temperature Heating Control of Air Flow in Diverging Part of Horizontal Nozzles, 
Indian Journal of Engineering, 2016, 13(34), 680-697, 


PageO 8 5 


ARTICLE 


4. NOZZLE SHAPE FORMULATION 


The nozzle inside shape is such that the area variation with distance along the nozzle is best given by [11] in linear curvature, we 


modified the area equation to be used for non-linear curvature with higher degrees as required. Simply we formulate the following 


equation, 


A(x) = a+ bx + cx? . (36) 


The values of coefficients a, b and c in equation (36) in terms of the inlet (throat) and exit areas will first be determined, 
When x=0.0 then A=A* 


When x=L then A=Ae 


When x=0.0 then dA/dx=0.0 


Applying these conditions into area variation equation, 


a=A*, at+bL+cL2 =Ae, b=0.0 


Equation (36) then becomes, 


. , (Ae — A*)x? 
A=A*+—* 7 .. (37) 
Now dividing equation (37) by A* to get the non-dimensional form of area variation, 
ER = 1+ (ERE — 1)X? (38) 


5. NUMERICAL FORMULATION FOR AREA VARIATION MODEL 


Because the equations governing fluid flow are non-linear, i.e. some of the terms in these equations involve products of flow 


variables; some form of iterative technique often has incorporated into solution procedure here [8]. 


Samheri, 


The flow domain is divided into equally spaced segments each of length Axwhere, 


Ax = —— (39) 


The conditions at throat section is specified, 
DR(1)=1.0, PR(1)=1.0, VR(1)=1.001, ER(1)=1.0 and ERE=ERE(N)=3.0 or 4.0. 


Near the throat (minimum cross-sectional area) the velocity ratio VR is so small that can be assumed as VR (i) =VR (1). This 
will make density ratio behave as, 


DR(1).ER(1 
DR(i) = ee Se ... (40) 


The results of density ratio are used from (2) to (N) points to find the values of pressure ratio of same points according to, 


PR(i) = PR(1 oy) : 41 
= PRO) (Seep (41) 
The velocities from (2) to (N)are updated using the first order finite-difference approximation to equation (35). 
The density ratio values are updated using equation (28) by employing under relaxation procedure to get, 
saGiaeue DRG)-VR@).ERQ) ie 
y= 1 r ~VR@.ER@ (i) roe ) 
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Then equations (11), (14), (23), (24) and (27) are used to find the remaining flow properties ratios. 


e = Acheck is achieved at all points for convergence of solution so that, 


IDR) new = DR(Doral 
DR@)oia 


< 0.001 (43) 


The estimation is carried out on Fortran Power Station using finite difference approximation to the space. In the present study air is 
chosen as the working fluid. 


6. NUMERICAL APPROXIMATION FOR ADIABATIC FRICTIONAL FLOW 


The effect of friction is usually introduced through momentum conservation [12], as shown in 


—dF — A.dP = p.v.A.dv (44) 
Frictional force will be found from, 
Tt. f.p. v2. A. dx 
—dF = SS a (45) 
VA 


Substituting equation (45) into equation (44) and dividing the result byp. Since dx represent the difference between any two nods its 
value is unity, 
T.f.v? dP 
ee Se ... (46) 
vA Pp 


Changing these values into dimensionless form, 


vi.f.v".VR? P*.dPR__ v*?.dVR? 


- (47 
‘A*.ER p*.DR 2 a 
Since, 
p* *2 
_—= * 2 = — — 
7 HR.T',M? = 1 =>, anda’ = ER(1) 
RT* = ve = a 48 
ae ..- (48) 
Equation (47) then becomes, 
T.f.VR2 dPR  dVR? 
a - = (49) 
ER(1).ER KDR 2 
Carrying out back spacing of the differentials, 
T.f.VRG)*  PRG)-PRG-1 
= te OF oko (i ) — vR@(VRG) —VRG- 1)) (50) 


ER(1). ER(i) k. DRG) 


Simplifying to the numerical model as, 


VR(i) = VRG—1) — _VTEVRG)_ _ PRG) — PRG — 1) (51) 


JER(1).ER@ ‘kK. DRG). VRG) 
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Then the effect of friction on the second term of right hand side in equation (51) and more exactly at average friction factor can 


be added. The values of average friction factor used in present study is chosen according to references [5], [6] and [9] as well as to 


experimental results from literature mentioned before in introduction. 


Most application values of f were restricted between 0.001-0.003 depending on surface roughness as well known. In present 


study we examine nine values of average friction factors. 


7. INTRODUCING OF HEAT TRANSFER AT ISOTHERMAL CONDITION FLOW INTOFRICTIONAL AREA 


VARIATIONFRICTIONAL MODEL FORMULATION 


The equation of state is, 


P=p.RT 
Its critical state is, 
P* = p".R.T* 
The dimensionless form is, 
PR= DR.TR 
Since for isothermal flow TR=1.0, 
PR=DR 
Dividing (i) state on initial state of equation (55), 
PR(i) = (GR) DR(i) 
DR(1) 
Since chocking takes place at (Zp velocity equation becomes, 
ne 
k kR.T* 
Simplifying, 
R.T* =v? 


The (k) values will be omitted from equations (49), (50) and (51) to be, 


vi.f£VRG)  PR(i) —PRG-1) 


VRG) = VRG—1)- = 


JER(1). ER(i) DRG). VR(i) 


The adiabatic law gives us the stagnation temperature ratio, 


The critical form will be, 


The dimensionless form becomes, 
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(57) 
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1+ (ae 
TZR = TR = (62) 
For isothermal condition, 
ta) 
T= a (63) 
VR(i) 
M= ... (64 
Vk (64) 
Entropy change equation becomes, 
S s=(<*)1 [=| RI E 65 
=\7-7) "| n| ps .. (65) 
In dimensionless form, 
s* =( <1 TR—InPR 66 
“hs =e n n .. (66) 
For isothermal flow, 
alee 67 
on n .- (67) 


8. ISOTHERMAL AREA VARIATION FLOW WITHOUT FRICTION MODEL FORMULATION 


The pressure ratio equation will be same as equation (56). Velocity ratio equation will be same as equation (34). Only Mach no. 
equation will be, 
M = VR(i) .. (68) 


Using equation (21), 


(69) 


The dimensionless Entropy change will be as equation (67). 


9. RESULTS AND DISCUSSION 


Extended results are found from the numerical analysis of mixed (isothermal-frictional-area variant) three types of flow with that of 
every double (isothermal-area variant) type or singular (area variant) type alone. Figures (2, 4, 6, 8, 10, 12, 14 and 16) show the flow 
properties ratios of pressure, velocity, density, area, stagnation temperature, Mach no., stagnation pressure and Entropy change as 
function of distance along nozzle length from throat till exit section for an area ratio configuration of 3.0. The difference is clear 
between each type of flow and a gradual change when the frictional forces are changed in the flow properties, which is too close to 
the experimental results from [15], as shown in figures (18-a &b).The same flow properties with Entropy change numerical results 
against nozzle length but with an area ratio configuration of 4.0 (for the same nozzle length as that of 3.0 area ratio means not just 
an extension) are shown in figures (3, 5, 7, 9, 11, 13, 15 and 17). A gain a clear change for the different flow types at different 
frictional forces. 

As well known the stagnation temperature and pressure as well as Entropy change are constant in non-frictional and non-heat 
transfer flow as shown in figures (10, 11, 14, 15, 16, and 17). 
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The pressure, density ratios and Entropy change in figures (2, 3, 6, 7, 16, and 17) show a clear decrease for all types of flow but in 


different amount. At the same time velocity, stagnation temperature and pressure ratios and Mach no. in figures (4, 5, 10, 11, 12, 13, 
14, and 15) show an increase in all types of flow and in different amount. 

A clear drop (difference) in pressure ratio (at exit section for its actual importance) exceed 10% between mixed and double flow 
type, and exceed 20% between mixed and singular flow type at 3.0 area ratio configuration, see figure (2). These differences at 4.0 
area ratio configuration dropped to 6% between mixed and double flow type, and dropped to 13% between mixed and singular flow 
type, see figure (3). Velocity ratio (at exit section) dropped to 70% between mixed and double flow type and 72% between mixed 
and singular flow type at 3.0 area ratio configuration figure (4). At 4.0 area ratio configuration the drop becomes 55% and 51%, 
respectively figure (5). Other differences in flow properties between the three flow types can be followed in different amounts as 
shown in the figures (6-17). 

In supersonic flow the increase of area will increase the velocity, while both the heat transfer (isothermal) and friction tend to 
decrease the velocity in different amounts. If we look at thedifferencesin properties throughout the nozzle, we will see that area 
variation was dominant at the first part of the nozzle on both isothermal and frictional forces. 

At the second part of the nozzle a reversal action takes place that isothermal and frictional forces become the dominant on area 
variation causing an inflection point either maximum or minimum depending on behavior of each property. Also it has been noticed 
that this inclination point is fixed (either minimum or maximum) for all properties at same frictional forces, as shown in table (1). 


Table 1 Isothermal-frictional-area variant points of inflection for all properties as function of nozzle length 


Area ratio f Minimum Point % Maximum Point % 
of Diverging Part of Diverging Part 
Length Length 
0.00100 91 70 
0.00125 89 66.3 
0.00150 87 62.6 
0.00175 85 58.8 
3.0 0.00200 83 55 
0.00225 81 51.3 
0.00250 79 47.5 
0.00275 77 43.8 
0.00300 75 40 
0.00100 95 75 
0.00125 93.2 71.3 
0.00150 91.3 67.5 
0.00175 89.4 63.8 
40 0.00200 87.5 60 
0.00225 85.7 56.3 
0.00250 83.8 52.5 
0.00275 81.9 48.8 
0.00300 80 45 
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Figure 4 Variation of velocity ratio along nozzle distance for different types of flow at area ratio of 3.0 
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Figure 7 Variation of density ratio along nozzle distance for different types of flow at area ratio of 4.0 
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Figure 8 Variation of area ratio along nozzle distance for different types of flow at area ratio of 3.0 
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Figure 9 Variation of area ratio along nozzle distance for different types of flow at area ratio of 4.0 
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Figure 10 Variation of stagnation temperature ratio along nozzle distance for different types of flow at area ratio of 3.0 
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Figure 13 Variation of Mach no. along nozzle distance for different types of flow at area ratio of 4.0 
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Figure 15 Variation of stagnation pressure ratio along nozzle distance for different types of flow at area ratio of 4.0 
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Figure 16 Variation of Entropy change along nozzle distance for different types of flow at area ratio of 3.0 
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Figure 18 (a) Experimental variation of static pressure ratio along nozzle distance, where 0.2-0.7 present subsonic back pressure 
ratios and 0.1present supersonic back pressure ratio [15], (b) Experimental variation of stagnation pressure ratio along nozzle 


distance [15]. 
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10. CONCLUSION 


The paper studied the full effect of friction at different degrees at the same time with heat transfer under isothermal condition in 


addition to area variation (mixed type flow) in non-linear curvatures. The bigger difference in flow properties is graduated from 
pressure, Entropy change, Mach no., stagnation temperature, velocity, stagnation pressure and density ratios, respectively. The 
numerical analysis gives a clear shift between the variables in the three types of flow (mixed, double and singular) and good 
agreement with the found experimental results. The new behavior of frictional-isothermal-area variant flow show a concave or 
convex shapes of variable changes, which are not known before either in double or singular flow types. 
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